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Abstract. We discuss the effect of curvature and matter inhomogeneities on the 
averaged scalar curvature of the present-day Universe. Motivated by studies of 
averaged inhomogcneous cosmologies, we contemplate on the question whether it is 
sensible to assume that curvature averages out on some scale of homogeneity, as 
implied by the standard concordance model of cosmology, or whether the averaged 
scalar curvature can be largely negative today, as required for an explanation of 
Dark Energy from inhomogeneities. We confront both conjectures with a detailed 
analysis of the kinematical backreaction term and estimate its strength for a multi- 
scale inhomogeneous matter and curvature distribution. Our main result is a formula 
for the spatially averaged scalar curvature involving quantities that are all measurable 
on regional (i.e. up to 100 Mpc) scales. We propose strategies to quantitatively evaluate 
the formula, and pinpoint the assumptions implied by the conjecture of a small or zero 
averaged curvature. We reach the conclusion that the standard concordance model 
needs fine-tuning in the sense of an assumed equipartition law for curvature in order 
to reconcile it with the estimated properties of the averaged physical space, whereas a 
negative averaged curvature is favoured, independent of the prior on the value of the 
cosmological constant. 



PACS numbers: 04.20.-q, 04.20.-Cv, 04.40.-b, 95.30.-k, 95.36,+x, 98.80.-Es, 98.80.-Jk 

1. The debate on the averaging problem and Dark Energy 

Homogeneous and isotropic solutions of Einstein's laws of gravitation do not account 
for inhomogeneities in the Universe. The question whether they do on average is a 
longstanding issue that is the subject of considerable debate especially in the recent 
literature ([Sj, [37] and follow-up references; comprehensive lists may be found in the 
reviews [27], [HE] and [H]). Averaging the scalar parts of Einstein's equations on space- 
like hypersurfaces of a foliation of spacetime pfl G3 [7] it was found that the Friedmannian 
framework is still applicable, however, one must include additional source terms due 
to the backreaction of inhomogeneities on a homogeneous-isotropic solution. These 
terms have geometrical origin and, as has been recently shown, can be represented by 
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a minimally coupled scalar field component, a so-called morphon field [15] . if those 
geometrical terms are interpreted as effective sources in a cosmological model with 
Friedmannian kinematics. This effective field can, like quintessence, other scalar field 
models [24], e.g. models motivated by higher-order Ricci curvature Lagrangians [25], 
[T6] or string-motivated effective actions [3], be employed to model Dark Energy. 
While the Newtonian and post-Newtonian frameworks suppress these effective scalar 
field degrees of freedom by construction [T3], [36], and so cannot lead to an explanation 
of Dark Energy, general relativity not only offers a wider range of possible cosmologies, 
since it is not constrained by the assumption of Euclidean or constant curvature 
geometry and small deviations thereof, but it is also needed to describe an effect 
that is strictly absent in a Newtonian model and in a standard (quasi-Newtonian) 
perturbation approach at a fixed background. This effect is reflected by the coupling of 
the fluctuations to the averaged model. In other words, fluctuations may be small, but 
measured relative to a non-Friedmannian background, and the evolution of this latter is 
most clearly expressed in terms of the evolution of effective geometrical properties such 
as the averaged scalar curvature we are considering in this paper (see [8] for detailed 
explanations). Speaking in favour of an averaged cosmology, it certainly enjoys the more 
physical status of incorporating inhomogeneities, and the clearcut fact that the effect of 
these inhomogeneities can be modelled by a scalar field speaks, by William of Ockham's 
razor, against introducing a dominating cosmological constant or an extra fundamental 
scalar field that is known to violate energy conditions in order to explain observational 
data. From this point of view one would also conclude that perturbation theory, 
if formulated at the background of a FLRW model with a dominating cosmological 
constant or an external scalar field source, would also not account for the physics behind 
the Dark Energy component. 

On the other hand, the FLRW cosmology provides a remarkably successful fitting 
model to a large number of observational data. As already mentioned, the price to pay 
is an unclear physical origin of either a dominating cosmological constant or an extra 
scalar field source that dominates recently. Given the fact that also a large amount of 
sources in the form of Dark Matter has not yet been detected in (non-gravitational) 
experiments, the standard model parametrizes an overwhelming fraction (95 percent) 
of physical ignorance. The generally held view, however, is that the FLRW cosmology 
indeed describes the physical Universe on average, which - if true - in turn asks for either 
a modification of the laws of gravitation, or the postulation of the above-mentioned dark 
sources of yet unknown origin. Moreover, the widespread use of the wording 'fitting 
model' is just name-dropping unless we devise a way to explicitly construct a smooth 
metric out of the inhomogeneous distributions of matter and curvature [261 [28] . In this 
- more refined - sense the FLRW cosmology is not a fitting model, rather it furnishes a 
conjecture on integral properties of the physical Universe that, as we believe, has to be 
first verified or falsified before more exotic vehicles of explanation are invoked. 

Both, the FLRW cosmology and a backreaction-driven averaged cosmology are 
candidates for the description of these integral properties, and in this paper we shall 
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estimate these properties from regionally (up to, say, 100 Mpc) observable quantities. 
For the FLRW cosmology in the form of the concordance model [39], [1], [611 [38] . 
the physical model is described by on average vanishing scalar curvature, while for 
a backreaction-driven cosmology, if we expect that Dark Energy can be fully routed 
back to inhomogeneities, the issue appears open. However, the consequences of a 
backreaction-driven model have been qualitatively, and to the extent we need also 
quantitatively, exploited in a number of recent papers (see [H] and references therein). 
For example, since a quantitative estimation of kinematical backreaction depends on 
specifying an evolution model for the inhomogeneities, the analysis of exact solutions 
like the Lemaitre-Tolman-Bondi solution (see, e.g., [2T], [5T], [22], [IS], and the reference 
lists in [56] and [8]), or scaling laws that satisfy the averaged equations [15] have been 
investigated. Consistency of an explanation of Dark Energy with the framework of the 
averaged equations has been demonstrated for both globally homogeneous cosmologies 
[15] (an assumption that we also adopt in the present paper) and, alternatively, 
globally inhomogeneous cosmologies [6], [7] . Although a quantitative evaluation of the 
backreaction effect in a generic inhomogeneous model is still to come, we already know 
a few features of a 'working' model, which are enough for our considerations. 




Figure 1. FLRW cosmology provides a successful fitting to observational data. The 
price to pay is an unclear origin of cither a dominant cosmological constant or of an 
extra scalar field source that appears to dominate the present dynamics of the Universe. 
Inhomogeneous hypersurfaces may be subjected to a smoothing procedure in order to 
find the corresponding smooth, i.e. constant curvature fitting model that we may call 
a FLRW template. We would then consider a hypersurface at a given instant of time; 
we cannot expect that the time-evolved inhomogeneous model can be mapped to a 
the constant-curvature evolution of the FLRW cosmology. 



We do not aim at investigating a 'fitting model' (see Figure 1) for the present- 
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day Universe in the strict sense mentioned above; this is the subject of ongoing work. 
We only remark here that any model for the evolution of the averaged variables can 
be subjected to a smoothing procedure in order to find the corresponding smooth, 
i.e. constant-curvature fitting model that we may call a FLRW template. We would 
then consider a hypersurface at a given instant of time; we cannot expect that the 
time-evolved inhomogeneous model can be mapped to a constant-curvature model. 
The above two candidates would provide different starting points, i.e. the initial 
data of a smoothing procedure for a given hypersurface are different. It was recently 
argued [52] that the averaged universe model (now both kinematically and geometrically 
averaged) could be represented by an effective FLRW metric with a time-scaling factor 
that differs from the usual global scale factor of a homogeneous-isotropic model and 
which is determined by the kinematically averaged Einstein equations. This ansatz 
for an effective metric assumes that smoothing (i.e. spatial rescaling) the actual 
matter and curvature inhomogeneities does not leave traces in the smoothed-out FLRW 
template metric at all times. In a forthcoming paper we are going to analyze this 
assumption in detail employing previous results on an explicit smoothing algorithm 
pTTj ITS] IT9] 13 QUI [HI dZ]. We emphasize that we are entitled to investigate integral 
properties of physical variables on a given hypersurface without entering the different 
question of whether this hypersurface (if actively deformed) can be effectively described 
by a 'best-fit' constant-curvature geometry. 

A rough guide that helps to understand the motivation of the present work is 
the following. For small inhomogeneities in the matter and curvature distributions, an 
approximate description of the cosmic evolution on average by a homogeneous solution 
of Einstein's laws of gravitation may be fine, but for the Late Universe featuring strong 
inhomogeneities, the validity of this approximation is not evident. We are going to 
address and justify this remark in the present work on the assumptions that (i) a 
homogeneous model satisfactorily describes the early stages of the matter-dominated 
epoch and (ii) there exists a scale of homogeneity. 

We proceed as follows. In Section [2] we look at the present-day Universe and 
device a three-scale model for it, where we hope that both readers, those who advocate 
the standard picture of the concordance model and those who advocate a backreaction- 
driven cosmology, agree. Then, in Section [3j we implement the details of this multi-scale 
picture and reduce the determining sources to those that are in principle measurable 
on regional (i.e. up to 100 Mpc) scales. Detailed estimates of the kinematical 
backreaction and averaged scalar curvature follow in Section [4j where we also provide 
simplified estimates in the form of robust bounds by, e.g., restricting the measurement 
of fluctuations to a comoving frame. In Section [5] we confront this latter result with 
the different assumptions on the actual averaged scalar curvature of the present-day 
Universe. 
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2. A fresh look at the present— day Universe 

2.1. Phenomenology: the volume-dominance of under dense regions 

Observing the Universe at low redshift returns the impression of large volumes that are 
almost devoid of any matter; a network of large-scale structure surrounds these voids 
that seem to be hierarchically nested and their sizes, depending on their definition, range 
from regional voids with less than 10 percent galaxy number content of the order of ten 
Megaparsecs [Mj . [30] . to relatively thinned-out regions of larger number density that, 
if smoothed, would span considerable volume fractions of currently available large- 
scale structure surveys. While the overall volume of the observable Universe seems 
to be dominated by underdense regions (the particular value of the volume-fraction 
of underdense regions being dependent on the threshold of this underdensity) , the 
small fraction of the volume hosting overdensities (groups, clusters and superclusters 
of galaxies) is itself sparsely populated by luminous matter and, this latter, appears as 
a highly nonlinear 'spiky' distribution. The phenomenological impression that matter 
apparently occupies a tiny fraction of space at all length scales could be questioned by 
saying that Dark Matter might be more smoothly distributed. Also, clusters contain 
a large amount of intergalactic gas, there are non-shining baryons, etc., so that the 
notion of an 'underdense region' has to be treated with care. However, simulations 
of Cold Dark Matter, assumed to rule the formation of large-scale structure, also 
demonstrate that voids dominate the present-day distribution. Again depending on 
the particular definition of a void, their fraction of volume occupation could, to give 
a rough value, be conservatively quoted as being 60 (23] percent in standard A— Cold 
Dark Matter simulations counting strong underdensities, and is certainly larger for more 
densly populated but still underdense regions. 

Thinking in terms of a homogeneous model of the Universe (not necessarily a 
homogeneous solution), i.e. a distribution of matter that on average does not depend 
on scale beyond a certain large scale (the scale of homogeneity), one would paint the 
picture of a redistribution of matter due to nonlinear gravitational instability. In a 
Newtonian simulation (where an eventually constant curvature of a FLRW spacetime 
is factored out on a periodic scale) this would happen in such a way that, due to the 
preservation of the overall material mass, an equipartition of overdense small-volume 
regions and underdense large-volume regions with respect to the mass content results, so 
that a sensible spatial average of the matter distribution must comply with the original 
value of the homogeneous density. In other words, the assumption that a volume- 
averaged distribution of matter would be compatible with a homogeneous model of the 
same average density seems to be a robust assumption, especially if inhomogeneities 
are dynamically generated out of an almost homogeneous distribution. This picture is 
true in Newtonian simulations, but for a subtle reason: although the time-evolution of 
the averaged density as a result of non-commutativity of evolution (time-scaling) and 
spatial averaging gives rise to kinematical backreaction, the periodic architecture and 
the Euclidean geometry of a Newtonian cosmology imply that these additional terms 
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have to vanish (see [13] for a detailed discussion of all these issues and proofs). In a 
general-relativistic framework this picture is in general false, even at one instant of time: 
the reason is that a Riemannian volume average incorporates the volume measure that 
is different for negatively and positively curved domains of averaging, and curvature 
does not obey an equipartition law (see Figure 2). 

Note also that a volume averaging on a Riemannian 3-surface could, even on the 
largest scale, introduce a volume effect in the comparison of the volume of a constant- 
curvature space and the actual volume of an inhomogeneous hypersurface (see [10J and 
[TT] for the definition and discussion of the volume effect; see also Hellaby's volume- 
matching example |32j). The standard model (but also the recent suggestion by [52J) 
implies that there is no such effect on large scales. It is illustrative to think of 2-surfaces, 
where curvature inhomogeneities always add up in the calculation of the total surface, 
so that there certainly is a large 2-volume effect due to surface roughening, but for 
three-dimensional manifolds, negative and positive curvature contribute with opposite 
signs, and so the 3-volume effect cannot easily be quantified. 




Figure 2. The picture that the scalar curvature in the physical space would average 
out on some large scale of homogeneity is naive in a number of ways. There is no 
equipartition law for the scalar curvature that would be dynamically preserved. 

Given the above remarks, one is no longer tempted to draw a picture of equipartition 
for the intrinsic curvature distribution. We shall, in this paper, not discuss the 
time-evolution of the scalar curvature (see [56] and [8] for detailed illustrations and 
discussions), large-time asymptotics [59j EQ], the role of a constant-curvature parameter 
in the fit to observations [22], or curvature models [501 EZ] (that are all related subjects 
of interest), but instead contemplate on the distribution of curvature at one given 
instant of time. Here, we demonstrate that the picture we would wish to establish 
in the concordance model, namely that the scalar curvature would average out on some 
large scale of homogeneity, is naive in a number of ways, and we shall implement the 
geometrical aspects of such a picture in Section [3j Obviously, this issue cannot be 
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addressed with Newtonian simulations; the curvature degree of freedom is simply absent. 
We know that in Riemannian geometry negatively curved regions have a volume that 
is larger than the corresponding volume in a Euclidean space section, and positively 
curved regions have a smaller volume, thus enhancing the actual volume fraction of 
under dense regions. 

We are now going to develop a multi-scale picture of the present-day Universe 
that is useful in the context of quantitative estimates, and that also helps to quantify 
multi-scale dynamical models, e.g. the one proposed by Wiltshire and his collaborators 
[62| [20| |63| I6H HOj [65], but it relates as well to any model that involves considerations 
of structures on different spatial scales. 



2.2. Scaling and coarse-graining: multi-scale picture of the Universe 

We are going to introduce three spatial scales. First, a scale of homogeneity Ln that 
could be identified with the size of a compact universe model, but need at least be larger 
than the largest observed typical structures. Second, a scale L £ that is as large as a 
typical void (within a range of values that depends on our definition of a devoid region), 
and third, a scale Lm that is large enough to host typical bound large-scale objects 
such as a rich cluster of galaxies. In observational cosmology we strictly have 

L n > L e > L M . (1) 

For the first we may also think of a length of the order of the Hubble-scale. The lower 
bound on this scale not only depends on the statistical measure with respect to which 
one considers the matter and curvature distributions as being homogeneous, but also on 
the concept of homogeneity that we have in mind. We here imply that averages of any 
variable beyond this scale will in practice no longer depend on scale, while generically, 
this may not happen at all. We do not claim here that this is indeed true, but we adopt 
this point of view in order to have a more transparent way of comparison with the 
standard model of cosmology. The assumption of existence of a scale of homogeneity 
may be a strong hypothesis; it is our choice of restricting the generality of the problem. 

According to what has been said above, we are entitled to assign different properties 
to the different scales, and we shall also sometimes idealize these properties in order to 
construct a simple but flexible model that reflects the phenomenology described above 
in terms of a small set of parameters. 

We start with an overview of the basic ingredients of our model and postpone 



details to Section [3j We employ the Hamiltonian constraint (see (13) below), spatially 
averaged on a given domain D (Scale Ly) that covers a union of underdense regions S 
(Scale Ls) and occupied overdense regions M. (Scale Lm)- We write the Hamiltonian 
constraint averaged over the first scale (for details see below and, e.g. [7]): 

(U) v = -ml -Q v + \QixG{q) v + 2A , (2) 

with the total restmass M-p := (q)v\T>\. The averaged spatial scalar curvature is denoted 
by {Ti)v, Hv '■= l/3(0)x> abbreviates the averaged rate of expansion 6 in terms of a 
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Figure 3. The three spatial scales to be used in our multi-scale picture of the Universe: 
the scale of homogeneity L-h, the scale of a typical void Lg, and the scale L_m hosting 
typical bound large-scale objects, e.g. a rich cluster of galaxies. 



volume Hubble rate, and the kinematical backreaction term Qx> encodes inhomogeneities 
in the extrinsic curvature distribution (or the kinematical variables); it is detailed in 
Section [3] (see Figures 3 and 4). 

Now, let us consider for illustrative purposes an idealization on the two other scales 
(in our concrete calculations later we shall indicate clearly when we make use of it): we 
require the volume Hubble expansion to be subdominant in matter-dominated regions 
and, on the other hand, the averaged density to be subdominant in devoid regions. In the 
first case, an expansion or contraction would contribute negatively to the averaged scalar 
curvature and so would, e.g., enhance a negative averaged curvature; in the second case, 
the presence of a low averaged density would contribute positively. We can therefore 
reasonably expect that, whether we use a strong idealization (see below) or a weaker 
distinction between over- and underdense regions, the overall argument based on the 
existence of such a partitioning enjoys some robustness. We shall also be able to condense 
our assumption on the partition between over- and underdense regions by introducing 
a parameter for the occupied volume fraction, X M := \V M \/\V\, where \T> M \ denotes 
the total volume of the union of occupied regions Ai; its value may be chosen more 
conservatively to weaken an eventually unrealistic idealization. At any rate, we shall 
keep our calculations as general as possible before we eventually invoke an idealization 
for illustrative purposes; it is only this latter quantity \m '■= \Dm\I\D\ that parametrizes 
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Figure 4. The fluctuations in the local expansion and shear generate a kinematical 
backreaction. This acts as a source term for scalar curvature when averaging the 
Hamiltonian constraint over the region of near-homogeneity T>. 



the (geometrical) volume-partitioning of the distributions of inhomogeneities. 

Note that, if we would strictly idealize voids to have (g)s = and and matter- 
dominated regions to have H M = for X M < A^, where A^ <C 1 is some critical 
scale (see Section [3] for more details on this transitory scale and the controling of this 
idealization), then we would have: 

(K)s = -m\ - Q E + 2A , (3) 

and 

(TZ) M = -Q M + 167iG(g) M + 2A , (4) 

together with 

H v = (1 - X M )H £ and (g) v = Xm(q)m ■ (5) 

This simplified view is useful as a rough guide on the sign of the averaged scalar 
curvature: consider for example the case where the kinematical backreaction terms 
in the above equations are quantitatively negligible, and let us put A = 0; we then infer 
that the averaged scalar curvature must be negative on Scale and positive on Scale 
L_A4, what obviously complies with what we expect. A non-vanishing A > is employed 
in the concordance model to compensate the negative curvature. 
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2. 3. Cosmological parameters and scalar curvature 

For our discussion we introduce a set of adimensional average characteristics that we 
define for the largest scale: 



OiJjy <Jiljy DlJj, Dflj, 

We shall, henceforth, call these characteristics 'parameters', but the reader should keep 
in mind that indexed variables are scale-dependent functionals. Expressed through 
these parameters the averaged Hamiltonian constraint (]2]) on the scale assumes the 
form: 

n£ + n? + + og = i . (7) 

In this set, the averaged scalar curvature parameter and the kinematical backreaction 
parameter are directly expressed through (TZ)x> and Qx>, respectively. In order to 
compare this pair of parameters with the 'Friedmannian constant-curvature parameter' 
that is the only curvature contribution in the standard model, we can alternatively 
introduce the pair (see, e.g., [7j) 

being related to the previous parameters by + &q n = fi^ + fig. 

For any of the smaller domains we discuss the corresponding adimensional 

parameters by dividing averages ( on the domain T always by H\ to avoid 

confusion. This will also avoid the pathological and useless definition of the cosmological 
parameters, e.g. on the domains M., where they are actually undefined in a strict 
idealization, since Hj^ is assumed to vanish. 

To give an illustration for the scale-dependence, note that, in the strictly idealized 
case, fi^ can be traced back to the average density in matter-dominated regions, 
(q)v — ^m(q)mi an d thus, inevitably, the density parameter constructed with an 
observed (q)m on the scale M. and divided by the global Hubble factor cannot be 
extrapolated to the global parameter. For example, a value today of a for this parameter 
would, for a volume fraction of matter-dominated regions of \m = 0.2, result in 

= 0.2a, i.e. a substantially smaller value that compensates the missing matter 
in the regions £, if they are idealized to be empty. Note in this context that a 
smaller mass density parameter on the global scale would also imply a smaller value 
of the necessary amount of backreaction. This can be seen by considering the volume 
deceleration parameter, 

^1=^ + 202-0? , (9) 

which, for A = 0, shows that decreasing the matter density parameter would also 
decrease the necessary backreaction in order to find, e.g. volume acceleration (q v < 0). 



Let us now discuss one of the motivations of the present work related to the Dark 
Energy debate. To this end, we have to explain why a substantial negative averaged 
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scalar curvature is needed, so that Dark Energy could be partly or completely routed 
back to inhomogeneities, i.e. fl^ = in the above equations, and, say, also « 0.2£j|] 
in conformity with the standard model of cosmology. First, by substantial we mean 
a cosmologically large negative curvature, i.e. expressed in terms of the cosmological 
parameters fi^° ~ 1. This parameter is exactly equal to 1, i.e. the averaged curvature 
is of the order of minus the square of the Hubble parameter, if there are no matter and 
no expansion and shear fluctuations. This applies to a simplified model of a void [56J. 
Second, in order to explain Dark Energy, the sum of the curvature and backreaction 
parameters has to add up to mimic a cosmological constant parameter of the order 
of « 0.7 in the case of a full compensation of the cosmological constant today, 

where we have adopted the assumption of the concordance model [39J of f2f = 0. A 
conservative scenario has been quantified in [TS] where it was assumed that the Universe 
at the Cosmic Microwave Background epoch is described by a weakly perturbed FLRW 
model and the amount of Early Dark Energy is negligible. The resulting scenario that 
would create enough Dark Energy features a strong curvature evolution from a negligible 
value to ~ 1, while the backreaction parameter must evolve from a negligible value 
to fig « —0.3, i.e. it has to be dominated by expansion fluctuations on the scale 
today. This scenario implies that the averaged scalar curvature must be close to the 
value of our simple void model, while at early times the averaged curvature parameter 
was compatible with zero. Speaking in terms of a morphon field [15], where the scalar 
curvature is associated with the potential of a scalar field, this scenario corresponds 
to the phantom quintessence sector (with negative kinetic energy of the scalar field), 
in which extrinsic curvature fluctuations (kinematical backreaction) grow slightly. We 
shall adopt this (present) value of the curvature parameter in our analysis as an extreme 
candidate for a backreaction-driven cosmology. 

Finally, we wish to emphasize that we are looking at spatial integral properties 
of inhomogeneous models of the Universe at the present time. This study does not 
include the important questions of (i) how the present-day structure we are looking at 
evolved dynamically out of an earlier state, and (ii) how these integral properties would 
relate to deep observations that necessarily involve considerations of the inhomogeneous 
lightcone and observable averages along the lightcone. We shall later propose strategies 
for the determination of these integral properties, also from observations. As a rule of 
thumb, observational results may be directly used in a shallow redshift interval, where 
the lightcone effect would be subdominant; galaxy catalogues as they are compared 
with the spatial distribution of fluctuations "at the present time" in simulations is an 
example). Our main result is a formula for the spatially averaged scalar curvature 
involving quantities that are all measurable on regional (i.e. up to 100 Mpc) scales, and 
it is therefore accessible by a shallow redshift interval. 

| We specify the index T>o as soon as we go to numerical estimates, indicating the value of the 
corresponding parameter today. Note that, in the averaged models, we may weaken this constraint, 
since we can allow for a scale dependence of this parameter in contrast to the situation in the standard 
model; see [5] and [33] for related discussions. 
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Eqs. (98) and (99), and their discussion thereafter. In particular, disregarding 



the contribution of gravitational radiation of cosmological origin, the adimensional 



backreaction term VLq can be written as (see (98)): 
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together with the formula for the adimensional averaged scalar curvature term: 



n v 



1 — f2™ — Vln with Vln from above 



VP 



V 



(10) 



where the rough meaning of the various terms is pictorially described in Figure 5. 




Figure 5. A pictorial representation of the meaning of the various terms characterizing 
the expression ( 10 ) for the adimensional backreaction fig. The Hubble parameters H-r 
and their variances S 2 H-r, as TZ varies in the various regions T>, Ai, and £ , are defined 
in § 3.3 The correlation lengths Lve,ve, Ljj, L^e,j and the parameters Lsh k , Vg[T>] 
are characterized in 5 14.21 



The following sections also prepare future work, e.g. on the determination of an 
"optimal frame" in which the variables are to be averaged in an evolving cosmological 
hypersurface. 
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3. Multi— scale analysis of curvature 

3.1. The Hamiltonian constraint 

In order to discuss the geometric structure of spatial curvature and of its fluctuations 
in observational cosmology, let us recall the essential steps required for constructing a 
cosmological spacetime out of the evolution of a Riemannian three-dimensional manifold 
E, which we assume for simplicity to be closed and without boundary. Note, however, 
that such a condition is not essential for our analysis and in due course it will be 
substantially relaxed. The geometry and the matter content of such a three-manifold 
is described by a suitable set of initial data (latin indices run through 1,2,3; we adopt 
the summation convention) 

(S , g ab , K ab , g, J a ) , (12) 

subjected to the energy (Hamiltonian) and momentum (Codazzi) constraints: 

1Z + K 2 — K\K b a = 16-kGq + 2 A ; V b K b a - W a K = 8nG J a , (13) 

where A is the cosmological constant, K := g ab K ab , and where 1Z is the scalar curvature 
of the Riemannian metric g ab ; the covariant spatial derivative with respect to g ab is 
denoted by V . 




Figure 6. The kinematics of the 3 + 1-splitting of a cosmological spacetime. The 
second fundamental form is here represented by the deformation of the shaded domain 
along the spacetime vector field t = N + N n, where n denotes the normal to E t , and 
N, N are the lapse function and and the shift vector field, respectively. 
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If such a set of admissible data is propagated according to the evolutive part 
of Einstein's equations (see Figure 6), then the symmetric tensor field K ab can be 
interpreted as the extrinsic curvature (or second fundamental form) of the embedding 
it : E — > of (£,<7 a &) in the spacetime (M^' ~ E x R,(jf( 4 ') resulting from the 

evolution of (T,, g ab , K ab , g, J a ), whereas g and J a are, respectively, identified with the 
mass density and the momentum density of the material self-gravitating sources on 
(£, gab). For short we shall call (E, g ab , K ab , g, J a ) the Physical Space associated with 
the Riemannian manifold (T,,g ab ) ■=— > (M^ ~ E x R, g^'). In what follows we shall 
make no use of the evolutive part of Einstein's equations and accordingly we do not 
explicitly write it down. 

3.2. Averaging scales 

Let us recall the various hierarchical length scales involved that we have described in 
Section |2j now associated with the curvature structure on the physical space (E,g a6 ): 
(i) the length scale L-^ defined by a spatial region V over which (E,g a b) can be 
viewed as describing to a good approximation a homogeneous and isotropic state 
(being not necessarily a homogeneous-isotropic solution of Einstein's equations); (ii) the 
length scales associated with the smaller domains over which the typical cosmological 
inhomogeneities regionally dominate, with an alternance of underdense regions and 
matter-dominated regions M (i \ with V = {U a £ (a) } U {U.^ (l) }, n M (i) = 0, 
£(«) n £ (/3) = 0, M {i) n M U) = 0, for all a ^ (3 and for all i ^ j. We denoted these 
latter length scales by and Lju respectively. For the former we sometimes say simply 
'voids' or 'empty regions', but our calculations are kept more general (see Figure 7). 



Underdense Regions 




Matter-Dominated Regions 



Figure 7. The alternance of Underdense and Matter-Dominated regions partitioning 
the spatial domain T>. 
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In order to discuss the implications generated by such a partitioning of £>, let us 
rewrite the Hamiltonian and the momentum constraints (13) over (T,,g a b) as 

K = 16ttGq + 2o 2 - -6 2 + 2A , (14) 

3 

Wo* = ^V b 9 - 8irGJ b , (15) 

where 6 := — K\ is the local rate of expansion, and a 2 := l/2a ab o- ab the square of the 
local rate of shear, with a ab := — (K ab — ^g ab K) being the shear tensor defined by 
-?T a b. We wish to average (14) over the region V = {U a £^} U {UjA^^} and discuss to 
what extent such an averaged constraint characterizes the sign of the curvature. Let us 
observe that on the scale of near-homogeneity the averaged Hamiltonian constraint 
(14) is assumed to have the structure Q 

(1Z) V = -6iZ| - Qv + 16itG(q)v + 2A , (16) 

where 

iZ2>:=^<0)j>, (17) 

is the (average) Hubble parameter on the scale L^. This assumption characterizes the 
kinematical backreaction term Qx> as 
2 

3~' 

Our strategy will be to express both (11) v and Qx> in terms of the typical local 
fluctuations of g, a, and 6 in the voids 8^ and in the matter dominated regions 
My i (see Figure 8). As a preliminary step, let us consider the average of a generic 
scalar-valued function / over V, 



Qv = ^(6 2 )v-QH^-2(a 2 ) v 



(fh:=\T>r / fdfi g , (19) 



v 



where \T>\ := j v dfi g , and dfi g := ^/gdx l dx 2 dx 3 , the Riemannian volume of T>. If we 
partition T> according to T> — {U a ^^} U {UiM.^}, where all individual regions are 
disjoint in the partitioning, then we can rewrite (f)v a s 

(f)v = \V\- ir £\E^\(f) eM + W-^IMVMmv . (20) 



where 



Since both the averages (f)s(.a), (f)_M^ an d the corresponding regions |£' a )| and |.M^| 
may fluctuate in value and size over the set of underdense {£^} a =i,2,... and overdense 
regions {M.^}i=i,2,...., it is useful to introduce the weighted averages of (/)g(a) an d °f 



: = ^ = i^r / , (2i) 
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\M^\(f) 



(f) 



M 



where \Vs\ 
\Ve\+ \T>m\ 



(f) 



EJ^ (a) l, and \V M \ 
T>\, we have 



YD 



M\ 



1-1 
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(22) 



23 



W\ 



J2i \M®\. Since V = {U a eW} U {U ?; .M«}, 

(23) 



If we now introduce the adimensional parameter 

\T> M \ 



A 



M 



W\ 



we can write (23) equivalently as 



(f)v = (1 - \m) (f)s + \ M (Dm- 



(24) 



(25) 




Figure 8. Splitting the average of the function / over the spatial region V into the 
weighted contributions coming from the underdense and from the matter-dominated 
regions. 



Applying in turn this formula to the volume-average of the scalar curvature 1Z, 
then we simply get 

(R) v = (1 - X M ) (K) e + Xm {K) m , (26) 



which, according to ( 14 ) implies 



1-A 



M, 



(n) v - 2A 



+ A 



M 



16nG( Q } M + 2(a 2 } M --(6 2 } M 
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3.3. The Kinematical Backreaction 

At this stage, we can look at the regional Hubble parameters H £ and H M assigned to 
the empty and matter-dominated regions and their associated mean-square fluctuations 
5 2 H £ and S 2 H M according to 

H £ := i H M := -(9) M ; 

S 2 H S := I {(6 2 } £ - (6)1) ; 5 2 H M := 1 {(6 2 } M - {Bf M ) . (27) 
One easily computes 

(K) v - 2A = (1 - X M ) [l^G(g) e - 6H 2 £ - Q £ ] + X M [16tiG(q) m - 6H 2 M - Q M ] , (28) 
where Q £ and Qm denote the kinematical backreaction terms on the respective scales: 
Q e := 6 5 2 H £ -2(a 2 ) £ (29) 

and 

Q M -=6S 2 H M -2(a 2 ) M . (30) 



If we insert (28) into the expression (16) characterizing Qx> we get 

Qv = {l~ ^m)Qs + \ M Qm + 6(1 -X M )H 2 + QX M H 2 M - 6H 2 . (31) 

Since 

H-p = (1 — Xm) H £ + Xm Hm , (32) 
a direct computation provides 

Qv = (1 - X M )Qs + X m Qm + 6Aa4(1 - A M ) (H e - H M f , (33) 
or, more explicitly, 

Q v = 6(1 - X M ) 6 2 H £ + QX M 5 2 H M + 6X M (1 - X M ) (H £ - H M f - 2 (a 2 ) v . (34) 

The above formulae for the averaged curvature and kinematical backreaction are general 
for our choice of a partitioning into overdense and underdense domains. 



It is important to observe that in the factorization (32) both H £ and Hm are 
effectively functions of A^, (this is simply a fact coming from the definition of the 
average factorization we have used), and for discussing the meaning of the expressions we 
obtained for the kinematical backreaction Q v and scalar curvature (TZ) x> it is often useful 
to assume a reasonable scaling for H £ (Xm) an d Hm(^m)- Our basic understanding is 
that, on small scales, the local dynamics of gravitationally bound matter will obliterate 
Hjn, whereas, if matter happens to be distributed over larger and larger domains, then 
it will more and more participate in the global averaged dynamics. By continuity, 
there should be a scale A^ marking a significant transition between these two regimes. 
Clearly, one can elaborate on the most appropriate model for such a transition, but the 
one described below, basically a Gaussian modeling, is quite general and has the merit 
of avoiding sudden jumps in the behavior of Hm- Also, it can be a natural starting point 
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for a more elaborate analysis. Thus, we wish to make an idealization by assuming the 
"stable-clustering hypethesis" to hold on the matter-dominated regions, H M = 0, it 
will only hold up to a critical scale X M G (0, 1), whereas for X M < Xm < 1 the quantity 
Hm(^m) smoothly increases up to H-p. To achieve this we model the scale-dependence 
of H M (X M ) according to 

(Aa< - l) 2 



H M (X 



M, 



Hx> exp 



(A- - 1)^ - (x M - iy 



for 



*M — 



< 1 



(35) 



and 



H M (X 



M, 



for < \ M < A 



(36) 

M £ [0, 1] with support 



It is easily verified that Hm(Xm) is a smooth (C°°) function of A 
in A^j < \m < 1) an d which vanishes together with all its derivatives at Xm = X M . (The 
graph of Hm(Xm) (see Figure 9) is the left-half of a bell-shaped function, smoothly 
rising from zero at Xm = X M , and reaching its maximum at Xm = 1). 



h m(^m) 




Figure 9. The graph of Hm(^m)- 



From the factorization Hx> = (1 — Xm) H £ + Xm H m we then get 



H £ (X M ) --=7^- l-A^e^te 5 ) for A^<A M <1, (37) 



H 

1 — Xm 

and 



H £ (X M ) := — \- for 0<X M <\ M • (38) 
i — am 

In light of the above remarks, we can, for < Xm — A^ , consider the idealization 

Hx> — (1 — Xm)Hs and Hm — 0, which together with (g)s = and (g)x> — Xm(q)m, 
allows to write the total kinematical backreaction in the simpler form: 

Qv = (1 - \ M )Qs + \mQm + -^—6Hl . (39) 

i — am 
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Thus, if we consider the fluctuations in the regional expansion rates 5 2 H^ and 5 2 Hm, 
the effective regional expansion rates themselves, as well as the regional matter contents 
Q £ m and and the volume fraction of matter X M , as in principle observationally 

determined quantities^ our main task will be to estimate the shear terms in VLq from 
its regional contributions in order to find an estimate for the curvature parameter. We 
shall now turn to this problem and investigate the shear terms that contain matter 
shear and geometrical shear and are therefore, like the intrinsic curvature, not directly 
accessible through observations; thus, we aim at deriving estimates that trace the shear 
terms back to the above more accessible parameters. 



§ Of course, the interpretation of these observations would still be model-dependent, sec the discussion 
in Sect. |5] 
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4. Estimating kinematical backreaction 

To provide control on Q-p we need to better understand the relative importance of the 
shear term (c 2 )v- Notoriously, this quantity is difficult to estimate without referring 
to particular models for the matter and gravitational radiation distributions, and in 
cases where we do not describe the variables in a comoving frame. In order to 
avoid a specific modelling and a specific coordinate choice, our strategy here is to 
express (cr 2 )v in terms of the correlation properties of the shear-generating currents. 
These correlation properties characterize the length scales over which the shear can 
be physically significant and allow for a natural parametrization of (o~ 2 )t>. Also, 
the underlying rationale is to separate, in each underdense and overdense region £ 
and Ai, the contribution to the shear coming from the local dynamics of matter 
and expansion, and to isolate this from the contribution coming from cosmological 
gravitational radiation. Technically, this procedure is based on a decomposition of the 
shear into its longitudinal (matter-expansion) and transverse (radiation) parts with 
respect to the L 2 (E) inner product defined by 

(W, V) L 2 (E) := J g ll g km W lk V lm d N , (40) 

where Wik, Vi m are (square-summable) symmetric bilinear forms compactly supported 
in E. With this remark along the way, we start by observing that for a given a ab we can 
write 



0~ab = 0-± ab + (X\\ ab , (41) 



with 



K^ll) L2(E) = 0, (42) 

where a± ab and a\\ ab respectively are the divergence-free part and the longitudinal part 
of a ab in E. 

Recall that a± ab generates non-trivial (proper time-) deformations of the conformal 
geometry associated with g ab , and represents the dynamical part of the shear, associated 
with the presence of gravitational radiation in E. Conversely, the term <j\\ a b is a 
gauge term associated with the infinitesimal Diff (E)-reparametrization of the conformal 
geometry generated by the motion of matter and inhomogeneities in the gradient of the 



expansion rate. The L 2 (E)— orthogonality (42) of these two terms implies that 



a dfj, g = / a± dfig + / cry d/j, g , (43) 
E Jt, Jt, 

where a ± 2 := l/2g ab g cd a ±ac (T ±bd and cry 2 := l/2g ab g cd a\\ ac op d . Explicitly, the terms a± ab 

and a\\ab (see Figures 10 and 11) are characterized by 

VV ±a6 = (44) 

2 

cr\\ab = V a w 6 + V b w a - -g a bV c w c ■= £wg a b , (45) 
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Figure 10. The longitudinal shear a\\ a i, associated with a 1-parameter family of 
volume-preserving diffeomorphisms in a fixed background geometry. The shearing 
diffcomorphisms are generated by a vector field w, solution of the elliptic equation 
(47 1, with sources given by §V&0 — 87rGJf,. 



where g ab denotes the conformal Lie derivative of the metric in the direction of the 
vector field w a . This latter is the solution (unique up to conformal Killing vectors-see 
below) of the elliptic partial differential equation V" (£wg a b) = VV a fc, i.e. 

Aw b + ^V b (V a w a ) + n ab w a = VV a6 . (46) 




Figure 11. The transverse shear a± a b associated with a volume-preserving 
(infinitesimal) deformation of the background metric geometry of E. 
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If we take into account the momentum constraints (15), V a cr a {, = ^VbO — 8irGJb, 



the elliptic equation determining the vector field Wb becomes 

A a b w a = ^V b e - SnGJb , (47) 

where 

A%:= A5 a b +^V a V b + n a b (48) 
is the formally self-adjoint elliptic vector Laplacian on the space of vector fields on 



4-1. The region of near-homogeneity T> 

In order to characterize geometrically the region of near-homogeneity D let us recall 
some basic properties of the vector Laplacian A\ introduced above. As x varies over 
(£,#), let {i^(i)(x)}i=i^,3 denote a local frame in T X S. For later convenience, we 
normalize the vector fields x \— > if(i)(x) according to J s ^)^a d/i g = and, 

when necessary, we assume that {ijj(i)(x)} 1=1^,3 is an orthonormal frame in T X H. For 
each given x 1— > ^^(x), let eu-\(x) be the corresponding vector field defined by the action 
of the operator A a b , i.e., 

C$ ef := A\ # , (49) 

where Cu\ is a normalization constant (with the dimension of a length squared) that 
will be chosen momentarily We denote by a® '■= . 9ab the shear associated with 

Observing that a^ b 'V a ip^ = \ a^ b a^, (no summation over (2)), and integrating 
by parts we get 

which, if we exploit the relation V a (a^ b ) = A a b ipa \ (see 46), yields 



2 dfJLg = 2 / (-A",) ^ dflg - (51) 



This latter shows that (— A a b ) is a positive operator whose kernel is generated by the 
conformal Killing vectors {£( a )} of (E,g a b). By proceeding similarly, from the relation 
° { 3 V% = \ . with i ^ k, we get 

J ^)V>2 } ) d» 9 = -\ jf a$ a* d N , (52) 



and 

''„/, u in "/'// ' ~ I ' 'in 
JIl, JY. 

In the absence of local conformal symmetries, which we assume to be generically the case 



< d »9 = 2 / 4) ^ d/ig ■ (53) 



in our setting, the left member of (51) is different from zero, and (49) can be interpreted 



as the statement that, if there are no conformal Killing vectors, then there always 
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exists a current e*(j) generating the given shear vector potential ip^y Moreover from 
V a (a { 2) = CT}e { a ] and |5o[ we get (dividing by Vol(Z,g)), 



(0\ _ 1 n. . I~. .2 



(i) e 6"7 E = -2 C '« ^ (<)2 ^ • (54) 

Note that (cr(j) 2 ) s is a natural measure of the characteristic length scale over which 
the conformal part of the metric g a b, (a tensor density of weight — |), changes along 
the vector field ipu\(x). One should also remark that the ratio 

(55) 



i i 

' 2 



[<<X(i) 2 >Ep [<*(*) 2 >e]" 

provides, for i ^ k, the typical correlation among the directional shears {o~^}k=i,2,3 
as the reference direction varies in {/0(i)}i=i,2,3- Even if the manifold (£,g) does not 
admit conformal symmetries, (since we assumed ker A b = 0), we should be able to 
geometrically select in (£,<?) a sufficiently large domain T> of near-homogeneity and 
near-isotropy. To this end, let us assume that in (X,g) there is a region (containing a 
sufficiently large metric ball) in which the local orthonormal frame x i— > {ip(i){x)} can 
be so chosen as to make (55) vanish, and the ratios 

ig °(i) 2 d ^g = ^ ) ( 56 ) 

h^ti)^ d ^9 

(where, in the latter expression, we have exploited the normalization of {^(i)}); are as 
small as possible with 



/ (J(l) 2 dflg = / (7(2) 2 dflg = / (T(3) 2 dflg . 

is Jt, 7s 



(57) 



In order to extract a natural geometric characterization of the region T> from these 
requirements let us recall that the first positive eigenvalue of the operator — A a b is 
provided by 

Al - mf l )' (58) 

where the inf is over all square-summable vector fields over (£,g) which are L 2 (X)- 
orthogonal to ker A a b . Since from (56) we have 

<^) 2 > s > Ai, 2 = 1,2,3, (59) 

this naturally suggests to select the local orthonormal frame {V'(i)} by choosing the vector 
fields as three independent eigenvectors, without nodal surfaces, of the operator 
— A a b associated with the first positive eigenvalue Ai, i.e., 

-A° 6 ^ i} = AxVfo, Jr^d^g = |E|5f , i,fc = l,2,3. (60) 

(we are here assuming that the multiplicity of Ai is > 3; for comparison, recall that the 
multiplicity of the first eigenvalue of the vector Laplacian on the standard 3-sphere is 6). 
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According to (49), with such a choice we have e^) = — AiV'(i) where the normalization 
constant takes the natural value CT^ := —Ax- These remarks allow to geometrically 

characterize the region of near-homogeneity as a region [V, {%jj^(x)} ) C (X,g), of size 



\V\ := A" 5 , (61) 

endowed with the minimal shear frame x i— > {ip®(x)} 1=1,2,3 generated by the chosen 
three independent nodal-free eigenvectors associated with Ai (see Figure 12). 




Figure 12. A pictorial representation of the (longitudinal) shear generated on a 
small sphere by three (independent) eigenvectors {4>(i)}i=i,2,3, associated with the 
first (degenerate) eigenvalue Ai of the vector Laplacian A^. Such eigenvectors define 
a moving frame in which the shear is, in a L 2 -average sense, minimized. The region 
of near-homogeneity D can be defined as the largest region in (E,g) where we can 
introduce such minimal shear eigenvectors. 



4-2. The shear correlation lengths 

We can infer a useful consequence of this geometrical characterization of T> if we 
exploit the Green function E%,(x,x'), of pole x G S, for the vector Laplacian A a b (see 
Figures 13 and 14). Let U C T> be a geodesically convex neighborhood containing 
the point x, (we can naturally restrict our attention to the region of near-homogeneity 
V C £). For any other point x' G U let l(x,x') denote the unique geodesic segment 
connecting x to x' . Parallel transport along l(x,x') allows to define a canonical 
isomorphism between the tangent space T X T> and T X >T>, which maps any given vector 
v(x) G T X V into a corresponding vector vp x , G T x /T>. If {C(h){ x )}h=i,2,3 G T X V and 
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{((k')(x')}k'=i,2,3 e T x iV, respectively, denote basis vectors, then the components of 
Vp iixx') can ^ e ex P resse d as 

k' 



,) (x f )=T k h '(x f ,x)v h (x), (62) 

where rfc (x f , x) denotes the bitensor e T x iT>®T*T> associated with the parallel transport 
along l(x, x'). This characterizes the Dirac bitensorial measure in U C T> according to 

b\ (x , x) := (x , x) S(x' , x) , (63) 

where S(x', x) is the standard Dirac distribution over the Riemannian manifold (T>, g) 
(see |11]). The Dirac bitensor 5^(x',x) so defined is the elementary disturbance 
generating the Green function for the vector Laplacian A a b according to 

x A\E^(x',x) = -6%(x',x), (64) 

where the pedex x denotes the variable on which the vector Laplacian is acting. In 
terms of E*(p, q) we can write 

^\x)=\ 1 [ E*(x',x)4)(x')d» g (x'), (65) 

which, since J" s ipa ^ (x) (x) d/J, g = |H| 8 % k1 yields 
Ai f f (fc) 



IE I 



^\x') E*'(x', x) r {l) {x) dii g (x') d N {x) = 6} . (66) 



This latter remark implies that 

(x' , x) d^ g (x') dfi g (x) := ^ ip%\x') E%'(x',x) ^(x) dfi g (x') dfi g (x) , (67) 

is a bivariate probability density describing the distribution of the geometric shear- 
anysotropies in (T>,g), along the minimal shear frame x i— > {ip(h)}h=i,2,3- If we define 
the components of the Green function with respect to {i>(h)}h=i,2,3 according to 

E*(x',x) := K'(x',x) r (i) {x) , (68) 



(the minimal shear Green function), then we can rewrite (67) in the simpler form 
p^(x', x) dfj, g (x') dfM g (x) := E^(x',x) d^ g {x') dfx g (x) . (69) 

It is with respect to the distribution of geometric shear described byp^)(x', x) dfi g (x') dfi g (x) 
that we can measure covariances and correlations of the shear generated, according to 
(47), by the distribution of currents §Vfc# — SnGJk- To this end let us consider the 
expression obtained from W^j, by contracting with the solution w b of (47) and inte- 
grating over (E,gab), i.e., 



J w b V a (a l]ab ) dp; g = ^J w b V b 9dfi g -8rrG J w b J b dfi g 



(70) 
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Integrating by parts, and exploiting again the identity cry a b V a w b = | a\\ a b cr^ b , we easily 
get 

2 



<7|| djJLg = 87rG W Jb djJLg — 



w° V b 9 dfj, 



9 ■ 



(71) 



(72) 



which, by taking the average with respect to (£,<?), provides 
<^|| 2 >s = ^G(w h Jbh-\{w h V b eh. 

In order to factorize the various contributions to (cr|| 2 )s coming from the various 
components of the current 8nG Jb — § V&6 1 , let us consider the integral representation of 
the solution of (47) along the minimal shear Green function (64), i.e., 

w a (x) = - / E^(x',x) j k (x') dfi g (x') 



(73) 



where, for notational ease, we have introduced the components jk{x') of the shear- 
generating current density with respect to the eigen-vectors basis {4>(k)}, according to 



j k (x') :=^ k) {x') [-V w 6-^kGJ w 



(74) 




Figure 13. The "current" j k (x') := 



<)(*') 



\Vh'9 — 8TrGJh<) x ,, at the point x' , 



generates a shearing vector field w(x), at the point x, by means of the Green function 
E^(x',x) of the vector Laplacian A^. In this illustration, the points x and x' are 
separated by a geodesic segment l(x,x'), and the shearing vector field w is shown as 
deforming a spherical domain, centered on x 7 into an ellipsoid. 



If we introduce (73) into (71) we get 

o\( L djjL g = / / j k (x') E^(x',x) j a (x)d/j, g (x')djjg(x) 



(75) 
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which, according to (|69|), implies 

jb(x') p^(af,x) j a (x) dfi g (x') dfi g (x) 



a\\ d\x g 



Note that the expression 



Gov (j, j) :-- 



Jb(x') pfa)( x 'i x ) f( x ) d/J,g(x) djJg(x) 



(76) 



(77) 



can be naturally interpreted as the covariance, in the minimal shear frame (T>, {i^u)}), 
between the values of the current j attained at the point x G T> versus the values 
attained in neighbouring points y ^ x, as these points vary over T>. Thus (76) can be 
rewritten as 

mi 



o~\\ 2 dfXg 



Ai 



Cov (j, j) 



or, more explicitly, in terms of the factors J and V#, (see (74)) 
IS 



cry dfJLg 



Ai 



-Cov(V9, V0) + (8tt Gf Cov(J, J) —Cov(V9, J) 

9 o 



(78) 



(79) 



where Cov(o,») is defined from (77) in an obvious way. 




Figure 14. The Green function E^(x',x) of the vector Laplacian A£ can be 
used to describe the curvature induced correlations among the distribution of the 
shear generating "currents" jk(x') and jk(x). By averaging these correlations over 
the domain of near-homogeneity one obtains the various correlation lengths, which 
describe the different contributions to the averaged longitudinal shear. 



As customary, we standardize (the absolute values of) these covariances between 
and 1 by dividing them by the appropriate powers of (| Jl 2 )^ and (| Wl 2 )^, (we are here 
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tacitly assuming that the main contribution to the shear generating currents comes from 
V, i.e., (|^| 2 ) E = (l^| 2 )x)' an d (l^| 2 )s = (l^^| 2 )©)- This defines the corresponding 
correlation coefficients 



Cov (V0, V0) 



Cov (J, J) 



O>rr(Y0.Y0):= /|t ^ |2 \ > ; Corr (J, J):- ^ ^ 



and 



Corr ye, J) :-- 



Cov ye, J) 



V 



(80) 



(wni (mi 

Note that < Corr (V0, V0) < 1, < Corr (J, J) < 1, whereas -1 < Corr (V0, J) < 

_ i 

1. The scale of T> is set by \ x 2 , (see (61), thus it follows that to such correlations we 
can associate the corresponding correlation lengths according to 



Corr ye, vey 


i 
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Corr (J, J) 


i 

2 


Ai 


; Lj tJ ■= 


Ai 





and 



Corr ye, J) 



which provide the length scales over which the distributions of V0, and J are significantly 
correlated and thus contribute appreciably to the matter-expansion shear cry 2 in T> (see 
Figure 15). 

In terms of these correlation lengths we can write, (dividing both members by |E| 
and assuming that (o"|| 2 ) E = (^H 2 )^) 



>p 9 
32ttG 



v 



(82) 



T 2 



we\ 2 )l <|j| 2 > 



V 



where =F is the sign of — Corr ye, J). To extract useful information from this relation 
let us remark that 



'j v (e-(e) v ) 2 d H 
j v ye\ 2 d N 



is the typical length scale associated with the spatial fluctuations of the expansion |0 
with respect to the average Hubble parameter in V. Since {{e 2 ) v — (6 1 ) 2 -,) = 9 5 2 H-p, we 
can write 

(ye\ 2 ) v = 9Ljl6 2 H v . (84) 

If necessary, this expression for (iVfll 2 )^ can be resolved into its finer components in 
the overdense and underdense domains into which T> is factorized. 
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Figure 15. A pictorial rendering of the various correlation lengths associated with 
the local distribution of the matter current J and the distribution of the gradient of 
the rate of expansion V0. 




' f e (9-(9) g f dfi 

fg \V6\ 2 dflg 



(85) 



which, as above, can be interpreted as the typical length scales associated with the 
spatial fluctuations of the expansion |# with respect to the local averages of the Hubble 
parameters in the overdense and underdense domains M. and £ (see Figure 16). Thus, 
if we factorize (IV^I 2 )^ as 

(\VB\X = X M J!™ 1 '}", {(*)„- «D' M ) 

+ ( 1 ~ A ^) Jr%2 (We-^l) > (86) 
we can express the term (IV^I 2 )^ as 

<|V^| 2 >^ = Lji M ((0 2 ) M -{0) 2 M ) + (1-\ M )LJI ({e 2 ) E -{9f e ) . (87) 
Since ((9 2 ) M - (6} 2 M ) = 95 2 H M and {{6 2 ) £ - {6%) = 9 5 2 H £ we can re-express Q as 
(| V 9\ 2 ) v = 9 LJ> V 5 2 H V = 9 \ M Lji M 5 2 H M + 9 (1 - X M ) Lj 2 £ 5 2 H £ , (88) 
which implies the useful relation 

b 2 H v = \ M 5 2 H M + (1 - \ M ) ^ S 2 H £ . (89) 

L SH M L 5H £ 
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Figure 16. A representation of the length scale L$h m associated with the spatial 
fluctuations of the local expansion |# with respect to the average Hubble parameter H-p 
in V, and its relation with the variance 8 2 Hxi of Hp and the gradient of the expansion 



V0. The resulting expression can be factorized, (see (88)), in its contributions coming 



from the underdense regions £ and the matter-dominated regions M. partitioning the 
domain of near-homogeneity T>. 



In order to deal similarly with the matter-current term } appearing in (82), 
it is profitable to rewrite it as 

(\J\ 2 ) V = *M (q)m + (1 - Xm) (Q)l , (90) 

\q)m \e)e 

where (g) M and (g) £ respectively denote the mass density in the overdense and 
underdense regions M. and £. Let us introduce the adimensional ratios, (we are using 
units in which c = 1), 



v*[£] ■= 



(b)m £* (er MW \M^\ 

(| j I 2 )aiW /\_n;f/„\2 ._ n\ „„,i ..2 ._(l J l 2 )fiW /._ n ;f/„\ 2 . 



Since V 2 [.M] and V^[S] are weighted averages of the (squared norms of the) typical 
velocity of matter in the overdense and underdense regions and £( Q ) we can 

naturally interpret V^ 2 [.M] and V^[£] as the (squared norm of the) typical velocity of 
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matter in the respective matter-dominated portion Ai and vacuum-dominated portion 
£ of V (see Figure 17). Thus, we write 



(\J\% = *m V*[M] (gf M + C 1 - x m) v;\£] (g) 



2 

S ' 



(92) 

with V?[.M] < 1, Vg[£] < 1. If we normalize this latter expression by (g) 2 ? then we can 
define the (squared) typical velocity of matter in the region T> according to 



V 



(q)t> 



\Q)v 



(i-\ M )vm 



{q)\ 



(93) 




Figure 17. The average velocities V^lyW] and V e [£] in the matter-dominated M. and 
in the underdense region £, respectively. These quantities are defined in terms of the 
corresponding average densities (g) M and (g) £ . 



Inserting these parametrizations into (82), and normalizing to the squared effective 
Hubble parameter H£, we eventually get 



H 2 



T 2 
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6 2 H 
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4-3. The transverse gravitational shear 

At this point, it is important to stress that the norm of the transverse part 2a ± 



2 



g ab g cd o-±acO~±bd, is not determined by (|47j). The term (a_i_ 2 )x> is associated with the 



presence of a non-trivial initial rate of variation of the conformal geometry of (T>,g), 
i.e., with initial data describing the presence of gravitational radiation of cosmological 
origin in (T>,g) (see Figure 18). We define the total energy density of gravitational 
waves in T> by 



oz7TLt 6 n v g crit 



3 H 2 

where g^ rit := g-^ is the density formally associated with the "critical density" of the 
standard Friedmannian model (in the region T> of near homogeneity). The ratio 

~TT~ := ^gw > (96) 

dcrit 

describing the relative strength of the energy density of gravitational waves with respect 
to the critical density, is the quantity conventionally used in cosmology for describing 
gravitational waves of cosmological origin. Thus, we parametrize the shear term (cr_i_ 2 )x> 
according to 

cf> - ( a± ^v (Q7) 
Ugw ~ 12H* ' { ' 




Figure 18. The transverse shear (a± 2 )v describes the rate of variation of the 
conformal geometry of (T>,g), and can be associated with the presence of gravitational 
radiation of cosmological origin in the region of near homogeneity (T>,g). 



Equipped with the above estimates, we now consider the backreaction term Q©, 
(divided by 6 if- 2 ,). From (34) we finally arrive at the key-result of this paper: 
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{H £ - H M ) 2 9T/2[ ^^w,j (S2nG (g) v L 8Hv (S 2 H V ) 
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j sh t 



H 2 



together with the formula for the averaged scalar curvature (see Figure 19): 

■Q 



VL^ + = 1 — Q% — f2o with from above . 



(98) 
(99) 




Figure 19. The kinematical backreaction, as described by (98 1, has, in the course of 
evolution, led to the emergence of a dynamical curvature term Vl^ (99 1 in the region 
of near-homogeneity (T>,g). 



4-4- Bound on kinematical backreaction in the matter-dominated Late Universe 

We are now going to invoke approximate assumptions in order to illustrate the result 
{(98, 99)}. At this stage the reader may critically compare our set of assumptions with 
the assumptions he would wish to make. At any rate, the following discussion employs 
simplifying assumptions and a more refined analysis has to take care of the neglected 
terms. 

At the stage we arrived at with the above formulae for the backreaction term and 
the averaged scalar curvature, it is clear what is the potential contribution of the various 
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terms involved. In particular, if, in line with the above analysis, we now assume that: 
(i) 



V e [£] = ; V e [M] = 



(100) 



namely that we are describing the region of near homogeneity T> in a frame comoving 
with mattei[jj] 

(ii) 



^GW « 1 



(101) 



i.e., absence of significant gravitational radiation of cosmological origin, then we finally 



obtain for the total kinematical backreaction parameter (|39|): 
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(102) 



This provides a reliable estimate in terms of the natural physical parameters involved. 
As expected, the shear terms (responsible for the factors proportional to L^q^q/ L 2 Hm 
and L^q ^q/ L 2 5He ) tend to attenuate the overall negative contribution of the fluctuations 
of the Hubble parameter. In a self-gravitating system with long-ranged interactions 
it is difficult to argue, if the V9-V9 current correlation length L(V0,V0) is smaller, 
larger or comparable to the expansion fluctuation length Lsh v - On the one hand, 
the above estimate indicates that the attenuation mechanism due to shear fluctuations 
cannot compensate for the global (i.e. on the homogeneity scale) negative contribution 
generated by 5 2 H£/H% } and 5 2 Hm/ H%>, since - due to the assumption of existence of 
a scale of homogeneity and, in addition, due to our implicit assumption of a globally 
almost isotropic state - the large-scale bulk flow must cease to display correlations and 
will be on some large scale subordered to the global Hubble flow pattern. This latter, 
however, could display large-scale fluctuations or not, and it is therefore to be expected 
that the ratio of the correlation lengths must be small, in conformity with our setup, 
only if there are significant differences between the globally averaged homogeneous state 
and a homogeneous solution which, this latter, features no fluctuations. On the other 
hand, we are confident that, on the void scale, the two correlation lengths are certainly 
comparable, and the shear term may also dominate over the expansion fluctuation term 
in Qs. This latter property was generically found in the Newtonian analysis [H] and 
it can be summarized by the expectation that the kinematical backreaction term would 
act as a Kinematical Dark Matter rather than as a Kinematical Dark Energy on the 
scale of voids (see [8] for a discussion), whereby on the global scale the domination of 
the expansion fluctuation term is possible and would then argue for an interpretation 
as Kinematical Dark Energy. 

|| These, henceforth neglected, terms would be relevant in the interpretation advanced by Wiltshire 
and his collaborators [SI HB [Ml [Ml HOI ES] . 
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It appears to us that this is a fixed point of the analysis, since in any case a non- 

L 2 L 2 

vanishing fraction X°' we and ^2' ve is going to add a negative contribution to Qt> (a 

positive contribution to VLq) and, therefore, neglecting this term will still allow us to 
provide bounds on the expected (or for some model prior necessary) fluctuations. As 
an estimate for the large-scale asymptotics we may, along the lines of this reasoning, 
propose the following simple formula for a rough, but to our opinion robust lower bound 
on a negative kinematical backreaction parameter, valid on the largest scales: 

- K < (1 - X M ) ^ + \ M gg* + A M (1 - \ M ) {Hs ~J m)2 .(103) 

With this formula we can bound a positive kinematical backreaction Qd from above, 
which provides information on the maximally expected backreaction and, in turn, on 
the maximally expected magnitude of the global averaged scalar curvature, to which we 
turn now. 



4-5. Bound on the averaged scalar curvature in the matter-dominated Late Universe 



Adopting the above restrictions of the general formula (98) we can write down an 



estimate for a global upper bound on the sum of a positive curvature parameter 
(corresponding to a negative averaged curvature) and the cosmological constant 
parameter, valid on the homogeneity scale. Using ([7]), (103), and fl^ = (1 — A^)f2^ + 
XM^m > we immediately get 
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The above formula demonstrates that the major players in a discussion on the maximal 
magnitude of the averaged scalar curvature will be (i) the volume fraction of occupied 
regions, (ii) the matter densities on the two regional scales, and (iii) the fluctuations 
in the volume-averaged Hubble parameters on the two regional scales, all normalized 
by the global volume Hubble rate that, this latter, can also be determined through 
measurement of regional parameters in view of H-p = (1 — \m)He + XmHm- AVe are 
now going to discuss this latter result more quantitatively. 



5. Concluding discussion 



First, an important remark that is relevant irrespective of whether we use a simple or a 
more refined estimating formula, and also irrespective of how exactly we determine the 
right-hand-side of, e.g. the inequality (104) from models or observations: in any case 
we would obtain a non-conclusive result for the value of the averaged scalar curvature 
itself due to the obvious degeneracy of the left-hand-side of Eq. (104). This latter 
depends on the value for the cosmological constant parameter that we have to choose ad 
hoc. For example, suppose that the right-hand-side of ( 104 ) (representing the measured 
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physical properties) would return a value 0.7 (now and in what follows concentrating on 
values today), we could not discriminate between the concordance model, i.e. adopting 

= 0.7, and at the same time = 0, and a backreaction-driven cosmology with an 
intrinsic physical curvature parameter =0.7 and vanishing cosmological constant. 
Such a value would therefore be compatible with two models of very different physical 
nature. An estimate would only be conclusive by setting a prior on the value of the 
cosmological constant. Only in the case where we would "exorcise" a cosmological 
constant, = 0, a measurement of the right-hand-side of (104) would be conclusive, 
since we could clearly discriminate between a zero-curvature universe model and a 
model with large curvature parameter. This degeneracy provides an obstacle especially 
in the expected cases of (i) the need of a dominating positive cosmological constant in the 
concordance model, and (ii) a substantial negative averaged curvature in a backreaction- 
driven cosmology (since in that case, by definition, > 0). However, despite this 
degeneracy, the physical status of the concordance model is weaker, since a conspiracy 
is required between the scale-independent value of the cosmological constant and the 
scale-dependent physical averages that can in principle be measured on any scale: even if 
such a conspiracy would hold on the largest scale T>, it is likely to be violated on smaller 
scales, a remark that essentially mirrors the status of a "fitting model" compared to a 
physical model for the averaged variables. 

Notwithstanding, an estimate of the right-hand-side of (104), or the more general 
formulae derived in this paper, is possible and useful, and we are now suggesting 
strategies for its realization. 



5.1. Strategy 1 - numerical simulations and analytical models 

Despite the fact that we wish to estimate Riemannian properties of the present- 
day Universe, the derived estimates for the large-scale averaged scalar curvature are 
determined by average properties on the scales of voids and matter-dominated regions. 
This in turn would suggest, as a first approximation, to determine the required regional 
parameters from a standard Newtonian N-body simulation. Thereby we accept to 
ignore the differences to the values that would be obtained for the Riemannian averages, 
but we expect to get an idea for the relative magnitudes of matter and expansion 
fluctuations within a well-studied framework. The determination of the regional 
parameters has to be done by controlled N-body simulations, this control respecting 
the mutual dependence of all the parameters involved: the volume fraction of matter- 
dominated regions, the calculation of averages and their fluctuation properties; both 
sensibly depend on spatial and particle resolution as well as on definitions of, e.g. 
void-finding algorithms and the employed thresholds on what we consider as over- 
and underdensities. In any case it is important to determine values within a single set 
of priors on the construction of the simulation and attaching a precise meaning to the 
involved parameters. 

Alternatively, an analytical estimate using relativistic models for the inhomo- 
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geneities can be investigated, first using perturbation theory along the lines of 
[371 El 1121 12], and its nonlinear extensions [53J, or studying constraints on the size 
of voids from swiss cheese-type or peak models [4~6l WA 158] . In this context we are 
currently generalizing the Newtonian analysis [H] , which is built on a relativistic gener- 
alization of a non-perturbative evolution model. In this model, the (root of the) metric 
fluctuations are calculated perturbatively, which is sensible, since the amplitude of met- 
ric perturbations is indeed small. Combining this with the exact averages then allows 
to take into account that the metric derivatives may not be small. 

In all these cases we can directly estimate the shear terms, since the assumption of 
a frame comoving with matter is adopted. 

Related to the time-evolution properties in N-body simulations or explicit 
analytical models we may also approach the problem dynamically from the point of view 
of effective (spatially averaged) Einstein equations. We may focus on the dynamical 
late-time properties of typical regions on the regional scales. As an example we 
may exploit the kinematically averaged Einstein equations jl], valid on any scale, and 
roughly estimate their time-asymptotics for a typical matter-dominated region. The 
aim would consist in understanding gravitationally bound systems and a corresponding 
virial equilibrium assumption that would constrain the involved energies. For example, 
if such a virialized state is characterized by a stationary "stable-clustering" volume Vm, 
then the averaged Raychaudhuri equation [lj on M. provides a relation between the 
averaged matter density and the kinematical fluctuation^} 

AnG(g) M = Q M + ^ M ; ^ = -4 (£1% + ^) . (105) 

Such a condition would replace the need for estimating fluctuations on typical domains 
Ai. With neglection of the shear term and terms contained in ^/j^ the above relation 
would read 4kG(q) m = QS 2 H M , or Q% = 245 2 H M , which still requires estimation of 

i 'm • 

5.2. Strategy 2 - observations 

The most direct, a priori model-independent strategy is, of course, to determine the 
values of the required parameters through observations. While this strategy is in 
principle free of priors on the spatial geometry, the interpretation of observational results 
has to involve model priors. Note that, at any rate, a model-dependence arises due to 
the fact that a volume-averaged value needs information on metrical properties, i.e. the 
volume depends on the geometry which is not directly observable (see, e.g., [331 H5] for 
strategies related to the determination of metrical properties). Actually all ingredients 
that are needed in our estimates are controversal in the literature, and accordingly an 

% We have inserted a phenomenological term together with the corresponding adimensional 

parameter, as an average over kinematical terms that include other forces relevant on regions A4. 
This term could have any sign and would be associated with velocity dispersion (compare |12j for a 
derivation of the corresponding non-averaged terms in a Newtonian setting - note also that vorticity 
is an important stabilization term that would have to be included in such a constraint). 
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observational determination would involve an ample range of values for the parameters 
and would suffer from the fact that different methods and interpretation mechanisms 
are necessarily involved and consequently would imply a loss of control on the mutual 
dependence of the parameters. 



5.3. Illustration 

Let us now put different model assumptions on the averaged curvature into perspective 
by just choosing some numbers for the physical properties on the regional domains. 



We write down the bound (104) for (i) the concordance model with (let us take rough 
numbers) = 0.3, = 0, = 0.7, and (ii) an extreme backreaction- driven 
model with = 0.3, = 1, and = 0; we obtain on the left-hand-side of jlQ4| 



the value LHS= 0.7 for the former and LHS= 1 for the latter. The right-hand-side of 



(104) reads: 



(H So - H Mo 



\2 



RHS = 1 - 0.3 + (1 - X Mo ) a e + ^m &m + \m (1 - Ax ) 9 - E ^ — , (106) 



with a £o := :i ^/^o) ; oc Mo : = (to) ■ 

Here, we have first ignored the possibility of choosing different values for the regional 
matter-density parameters to ease the comparison. To make the average expansion 



properties concrete, we employ our scaling models for Hs and Hm , see Eq. (35)ff. 



with a transition value XJ 4o =0.1. 

The reader may now specify the free parameters involved (A» , ae , a^J. 
Considering a maximal range for the volume fraction, < \m < 0-5, we then 
obtain RHS = 0.7 + a>s for a completely empty universe model, and RHS = 0.7 + 
0.5(a,f o +o>M ) +0.25 • 0.09, for a volume-equipartition of over- and underdense regions. 
This shows how global constraints imposed by our model priors would constrain the 
fluctuation terms. E.g., a value of ag = 0.3 would be the minimally necessary 
fluctuation amplitude of the expansion rate on the void scale, if the expectation from 
the extreme backreaction-driven model, that has practically emptied the Universe, is 
met. For a non-zero A» we have to also specify a_M - 

As another example let us now assume different density parameters on the different 
regions. For a range of values for the volume fraction of occupied regions Xm = 
0.2; 0.3; 0.4 (remember that this parameter is the relativistic volume fraction), we for 
example get with the priors on the global density parameter, w 0.3, and the density 
in voids = 0.03 (i.e. roughly 10 percent of the average density is found in voids) the 
corresponding values = 1.38; 0.93; 0.71. By working out a stationarity assumption 



in line with (105), we could then find a_M directly from the density parameter in 



matter-dominated regions. 

At this stage we leave it to the reader to exploit the result further and to understand 
the relations of assumed or measured numbers to the global model priors according to 
his/her experience. 
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5.4- Concluding remarks on the concordance model 

In order to reconcile the standard concordance model with the estimated bounds on the 
averaged scalar curvature, the above discussion indicates that this would be difficult, 
since for an assumed cosmological constant parameter of 0.7, the curvature parameter 
acquires additive pieces due to the fluctuations of the Hubble parameter on the two 
regional scales. Overall, the above (rough) discussion would suggest that, independent 
of the prior on a specific value of the cosmological constant parameter in between 
and a value that compensates the matter density parameter (in our example 0.7), the 
averaged curvature would be negative. 

The remaining question is, how robust this latter result is, given our assumptions. 
For this purpose we are taking the role of advocating the standard concordance model 
by pushing the estimates to their extremes. For this end we have to argue that all the 
additional contributions due to inhomogeneities are negligible, so that our formulae in 
fact reduce to the standard Hamiltonian constraint for a homogeneous distribution of 
matter and curvature. Only in this case the concordance model, i.e. a zero-curvature 
universe model, can be reconciled with the physical space. This statement implies that 
there must be a conspiracy between (Al) the shear terms on both regional scales, (A2) 
the expansion terms on both regional scales, and (A3) the difference between the Hubble 
rates on those regional scales, such that the overall contribution of these terms cancels 
on the homogeneity scale. Alternatively, we could conjecture that (Bl) the magnitudes 
of the involved fluctuation amplitudes in our formulae are quantitatively negligible and 
(B2) that Hs = Hm- However, these latter options are unrealistic, since we would have 
to require negligible fluctuation amplitudes also on regional scales (ad Bl), and that 
e.g. a cluster of galaxies participates in the full Hubble flow (ad B2). 

Whether the above conspiracy (A) could hold is difficult to test. However, we know 
that it holds in a Newtonian model due to boundary conditions [13], [TJ], i.e. this 
conspiracy is suggesting that the present-day Universe can be effectively described in 
Newtonian terms [US], i.e. it is equivalent to saying that the structure formation history 
had no impact on the evolution of the averaged curvature, and fluctuations have been 
decoupled from the geometry (being a dynamical variable in general relativity) until 
today. Although we cannot exclude this possibility, it relies on a fine-tuning assumption 
that can be physically justified for a Euclidean universe model, where the curvature 
distribution trivially satisfies an equipartition law (since there is no curvature). As we 
have discussed, we cannot expect such an equipartition of curvature in the general case, 
especially when the partitioning is biased towards a volume-dominance of underdense 
regions (negative curvature). 

Thus, to reconcile the standard concordance model with the physical average of the 
present-day Universe is essentially equivalent with postulating an equipartition law for 
the scalar curvature on the homogeneity scale T>, hence a fine-tuning assumption on the 
curvature distribution. Here, one should be aware of the fact that such an equipartition 
law must be the result of a conservation law for the scalar curvature in a dynamical 
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situation, similar as an equipartition of the density distribution on V being the result of 
restmass conservation. We know that such a conservation law does not exist in general. 
In the case of a dust matter model it is replaced by the condition (see [8] for a detailed 
discussion): 



i.e. a particular dynamical combination of kinematical fluctuations and averaged scalar 
curvature is conserved, not the averaged scalar curvature itself, that here would be 
represented by the conservation of the Yamabe functional (1Z)t)V^ 3 . This does not mean 
that such a conservation law is excluded, but it would be equivalent to an uncoupled 
evolution of fluctuations in a constant-curvature space section, since d t ( (TZ)T>V^ 3 ^j = 

implies (lZ)v V v 2 ^ 3 . It is interesting that the mechanism of a backreaction-driven 
model relies on a coupled evolution between Q v and (7Z)t> (thoroughly discussed in 
[3 US])) hence, it genuinely violates an equipartition of curvature. 

5.5. Summary 

We conclude that we need a fine-tuning assumption on the scalar curvature distribution 
(an equipartition law on the scale of homogeneity) in order to reconcile the concordance 
model as a viable model for the physical properties of the present-day Universe. 
Dynamically, this hypothesis implies an uncoupled evolution of kinematical fluctuations 
and intrinsic curvature. If such a hypothesis is not adopted, our investigations point 
to an overall negative averaged scalar curvature. Measuring fluctuations in the volume 
Hubble rate on the two regional scales together with the volume fraction of occupied 
regions would allow us to support or to rule out a large negative averaged scalar 
curvature, required for an extreme backreaction-driven model, only if A = 0. The crux 
in this consideration is the degeneracy by assuming a non-zero cosmological constant. 
However, if the large value for the fluctuations needed for a backreaction-driven model is 
achievable (for calculations and discussions of estimates see [H]; [56], Sect. 3. 3; (53] and 
|42| |4"3]). there would be good reasons to shift our interpretation from the assumption 
of a "curvature-compensating" A to the acceptance of a substantial negative physical 
curvature. 
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